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Relativistic many-body perturbation theory is applied to study properties of ions of the francium 
isoelectronic sequence. Specifically, energies of the 7s, 7p, 6d, and 5/ states of Fr-like ions with 
nuclear charges Z = 87 — 100 are calculated through third order; reduced matrix elements, oscillator 
strengths, transition rates, and lifetimes are determined for 7s — 7p, 7p — 6ci, and 6d — 5/ electric- 
dipole transitions; and 7s — 6d, 7s — 5/, and 6/5/2 — 6/7/2 multipole matrix elements are evaluated 
to obtain the lifetimes of low-lying excited states. Moreover, for the ions Z = 87 — 92 calculations 
are also carried out using the relativistic all-order single-double method, in which single and double 
excitations of Dirac-Fock wave functions are included to all orders in perturbation theory. With 
the aid of the SD wave functions, we obtain accurate values of energies, transition rates, oscillator 
strengths, and the lifetimes of these six ions. Ground state scalar polarizabilities in Fr I, Ra II, 
Ac III, and Th IV are calculated using relativistic third-order and all-order methods. Ground state 
scalar polarizabilities for other Fr-like ions are calculated using a relativistic second-order method. 
These calculations provide a theoretical benchmark for comparison with experiment and theory. 

PACS numbers: 31.15.Ar, 31.15.Md, 32.10.Fn, 32.70.Cs 



I. INTRODUCTION 

A detailed investigation of radiative parameters for 
electric dipole (El) transitions in Fr-like ions with Z = 
89-92 was presented recently by Biemont et al. jl|. The 
electronic structure of Fr-like ions consists of a single nl 
electron outside of a core with completely filled n—1, 2, 
3, 4 shells and 5s, bp, 5d, 6s, and 6p subshells. In Fig. [U 
we plot one-electron DF energies of valence 5/, 6d, and 
7s states as functions of Z. We find that the valence 7s 
orbital is more tightly bound than the 5/ and 6c? orbitals 
at low stages of ionization {Z = 87-89), while the 5/ 
and 6d orbitals are more tightly bound for highly ionized 
cases {Z > 90). Competition between the 5/, 6d, and 
7s orbitals leads to problems for calculations, making it 
difhcult to obtain very accurate excitation energies and 
line strengths for the transitions between the low-lying 
5/, 6d, and 7s states. 

Relativistic Hartree-Fock and Dirac-Fock atomic struc- 
ture codes were used in Ref. [l| to perform calculations 
of radiative transition rates and oscillator strengths for 
a limited number of transitions using the energies given 
by Blaise and Wyart [2| , where experimental values were 
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given for the 7p, 8p, nd (with n =10-33), and ns (with 
n =12-31) levels of neutral Fr, for 24 levels for Fr-like 
Th, and for seven levels of Fr-like Ac and U. Adopted 
energy level values in the ns, np, and nd (n < 30) se- 
ries of neutral francium were presented by Biemont et al. 
Q . Experimental measurements of energy levels of the 8s 
and 7c? states in Fr I were reported recently in Refs. 0, . 
The experimental energies of 13 levels of Fr-like Ra were 
reported in the NIST compilation 6]. 

Lifetime measurements for neutral francium were pre- 
sented in Refs. 0, i, i, [M [H for the 7p, 6c?, 9s, and 8s 
levels. In those papers, experimental measurements were 
compared with ab-initio calculations performed by John- 
son et al. 12], Dzuba et al. 13,[l3|, Safronova et al. [Tsj . 
and Safronova and Johnson 16| . Third-order many-body 
perturbation theory was used in Ref. [l^ to obtain the 
El transition amplitude for neutral alkali-metal atoms. 
The correlation potential method and the Feynman di- 
agram technique were used in Refs. [H, [l3| to calculate 
El matrix elements in neutral francium and Fr-like ra- 
dium. Atomic properties of Th IV ion were studied by 
Safronova et al. [17] using a relativistic all-order method. 

The present third-order calculations of excitation en- 
ergies of the 7s, 7p, 6c?, and 5/ states in Fr-like ions with 
nuclear charges Z = 87 — 100 start from a closed-shell 
Dirac-Fock potential for the 86 electron radon-like core. 
We note that Th IV is the first ion in the francium isoelec- 
tronic sequence with a [Rn] 5/5/2 ground state instead of 
the [Rn]7si/2 ground state as for Fr I, Ra II, and Ac III. 
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FIG. 1: Dirac-Fock one-electron energies {E/(Z — 84)^ in 
cm~^) for the 5/j, 6dj, and 7sx/2 states of Fr-like ions as 
functions of Z. 



Correlation corrections become very large for such sys- 
tems as was demonstrated by Savukov et al. 18], where 
the ratio of the second-order and DF removal energies for 
the [Xe]4/5/2 ground state in Co IV and Pr V were 18% 
and 11%, respectively. 

In the present paper, dipole matrix elements are cal- 
culated using both rclativistic many-body perturbation 
theory, complete through third-order, and the rclativis- 
tic all-order method restricted to single and double (SD) 
excitations. Such calculations permit one to investigate 
the convergence of perturbation theory and estimate the 
theoretical error in predicted data. To obtain lifetime 
predictions, multipole matrix elements for 7s— Qd^ 7s — 5 f 
and 5/5/2 ~ 5/7/2 transitions are also evaluated. Addi- 
tionally, scalar polarizabilities for the 7si/2 ground state 
in Fr I, Ra II, and Ac III are calculated using rclativistic 
third-order and SD methods. Finally, scalar polarizabili- 
ties of the 5/5/2 ground state of Fr-like ions with nuclear 
charge Z = 91-100 are calculated in second-order MBPT. 



II. THIRD-ORDER AND ALL-ORDER MBPT 
CALCULATIONS OF ENERGIES 

We start from the "no-pair" Hamiltonian [l^ 

H = Ho + Vi, (1) 

where Hq and Vj can be written in a second-quantized 
form as 



Ho = Sialai , 



(2) 



ijkl 



(3) 



Negative-energy (positron) states are excluded from the 
sums. The quantities Si are eigenvalues of the one- 
electron Dirac-Fock equations with a frozen core and 
gijki is a two-particle Coulomb matrix element. Our cal- 
culations start from a V'^^^ DF potential for a closed- 
subshell radon-like ion. 

The all-order single-double (SD) rnethod was discussed 
previously in Refs. 0, [13, [SI, [H, [2II2II . Briefly, we rep- 
resent the wave function of an atom with one valence 
electron atom as "ify = "^fP with 



SD 



E 



mnab 



+ E 



Et t 

m7ia 



^v, (4) 



where is the lowest-order atomic wave function, which 
is taken to be the frozen-core DF wave function of a state 
V. Substituting the wave function ^E'^'^ into the many- 
body Schrodinger equation, with Hamiltonian given by 
the Eqs. (IT]-[2]), one obtains the coupled equations for 
the single- and double-excitation coefficients Pmv, Pma, 
Pmnva, and Pmnab- Thc couplcd cquatious for the excita- 
tion coefficients are solved iteratively. We use the result- 
ing excitation coefficients to evaluate multipole matrix 
elements and hyperfine constants. This method includes 
contribution of important classes of MBPT corrections 
to all orders. 

The SD valence E^^ energy does not include a all 
third-order MBPT corrections. The missing part of the 
third-order contribution, is written out in Ref. [2^ 

and must be calculated separately. We use our third- 

(3) 

order energy code to separate out i?cxtia add it to 
the . For notational simplicity, we drop the index v 
in the designations in the text and tables below. 

Results of our energy calculations for low-lying states 
of Fr I-U VI are summarized in Table HI Columns 2-7 of 
Table [T] give the lowest-order DF energies -B^^-' , second- 
and third-order Coulomb correlation energies, and 
E'^^\ first-order Breit contribution B'^^\ second-order 
Coulomb-Breit B^^^ corrections, and the Lamb shift con- 
tribution, i?LS- The sum of these six contributions is 
our final third-order MBPT result E^^^^ listed in the 
eighth column. First-order Breit energies (column B^^' 
of Table m include retardation, whereas the second-order 
Coulomb-Breit energies (column B^^' of Table|I| are eval- 
uated using the unretarded Breit operator. We list all- 
order SD energies in the column labelled E^^ and the 
part of the third-order energies omitted in the SD cal- 



culation in column E, 



(3) 



We note that includes 



E^'^\ part of E''^\ and dominant higher-order correc- 
tions. The sum of the six terms E'^^, E^^^La: 
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TABLE I: Zeroth-order (DF), second-, and third-order Coulomb correlation energies E^"\ single-double Coulomb energies E^^, 
El'^^J.^, first-order Broit and second-order Coulomb-Brcit corrections _B'"' to the energies of Fr-like systems. The total energies 
are = i?*"' + E^''^ + + B^'^ + B*^' + Els, E^° = + + 75(3)^^ + ^d) + 5(2) + Els- Units: cm-^. 



nlj 


£;(o) 


£;(2) 


^(3) 


5(1) 


5(2) 




Els 


^ tot 




extra 


TT'SD 

^ tot 












Fr I, Z = 


= 87 












7Sl/2 


-28767 


-4763 


1737 


67 


-131 




13 


-31845 


-4658 


724 


-32753 


7351/2 


-18855 


-1847 


546 


30 


-35 







-20162 


-1991 


255 


-20597 


7p3/2 


-17655 


-1346 


391 


19 


-30 







-18621 


-1433 


184 


-18915 


64/2 


-13807 


-2424 


711 


19 


-62 







-15563 


-3180 


326 


-16990 


64/2 


-13924 


-2240 


616 


15 


-61 







-15594 


-2854 


288 


-16744 


8Sl/2 


-12282 


-1053 


396 


17 


-32 
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-12952 


-925 


162 


-13058 


8pi/2 


-9240 


-554 


171 


10 


-13 







-9625 


-551 


77 


-9716 


8p3/2 


-8811 


-421 


128 


7 


-11 







-9108 


-421 


58 


-9178 


74/2 


-7724 


-981 


298 


10 


-30 







-8427 


-991 


132 


-8604 


74/2 


-7747 


-860 


244 


7 


-28 







-8384 


-856 


109 


-8515 
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= 88 












7si/2 


-75898 


-7529 


2896 


147 


-250 




33 


-80634 


-6692 


1152 


-81508 


64/2 


-62356 


-8727 


2764 


155 


-398 







-68562 


-8042 


1152 


-69488 


64/2 


-61592 


-7537 


2202 


114 


-360 







-67174 


-7034 


926 


-67947 


7pi/2 


-56878 


-4182 


1370 


102 


-109 







-59698 


-4027 


587 


-60326 


7p3/2 


-52906 


-3130 


1011 


63 


-90 







-55053 


-3020 


433 


-55519 


5/5/2 


-28660 


-2563 


824 


11 


-63 







-30452 


-4438 


371 


-32780 


5/7/2 


-28705 


-2491 


784 


8 


-61 







-30466 


-4159 


353 


-32564 
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7Sl/2 


-133640 


-9552 


3739 


233 


-357 




58 


-139519 


-8192 


1456 
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64/2 


-130697 


-11506 
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-138927 


-10036 


1479 


-139617 
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-2747 
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1663 
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-458 
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-272 
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-153572 
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184 


-226 
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-5619 
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-336671 


-34071 
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-3390 







-360855 


-27589 
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-361865 


5/7/2 


-331505 


-32627 
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-26571 
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-297300 
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3556 


438 
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1446 


-307768 


7Si/2 


-274949 


-12646 


3359 


425 


-557 




126 


-284242 


-10606 


1348 


-284212 


7pi/2 


-232148 


-9630 


2566 


417 


-356 







-239151 


-8533 


1040 


-239580 


7j>3/2 


-216015 


-7600 


4939 


258 


-297 




1 


-218715 


-6925 


1900 


-221079 
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5/5/2 


-481613 


-35163 


12554 


1194 


-3929 







-506958 


-28393 


4876 


-507866 


5/7/2 


-474700 


-33781 


11863 
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-3789 







-499509 


-27398 


4604 


-500385 


64/2 


-404911 


-15830 


4675 


742 


-1275 







-416599 


-13737 


1908 


-417273 


64/2 


-396467 


-14000 


3718 


557 


-1181 







-407374 


-12472 


1538 


-408025 


7Sl/2 


-357141 


-13966 


5366 
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-655 




173 
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-11752 


2079 


-366763 


7pi/2 


-306870 
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3742 


548 
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FIG. 2: Coulomb second and third-order energies as functions of Z for the 5/j, 6dj, 7si/2, and 7pi/2 states in Fr-Uke ions. 



and Ei^s gives the final all-order results Ef^ listed 
in the eleventh column of the table. 

As expected, the largest correlation contribution to 
the valence energy comes from the second-order term, 
E'^'^\ This term is simple to calculate in comparison 
with E^^^ and terms. Thus, we calculate the E^"^") 
term with higher numerical accuracy than £'(3) and E^^ . 
The second-order energy E^"^^ includes partial waves up 
to Imax — 8 and is extrapolated to account for con- 
tributions from higher partial waves (see, for example, 
Refs. [1^, d^). As an example of the convergence of 
E^"^^ with the number of partial waves we consider the 
5/5/2 state in U VI. Calculations of i?'^) with Imax — 
6 and 8 yield i^'^' (5/5/2) = -33598 and -347559 cm-\ 
respectively. Extrapolation of these calculations yields 
-35163 and -35227 cm^^, respectively. Therefore, we es- 
timate the numerical uncertainty of £'^'^■'(5/5/2) to be 
64 cm~^. This is the largest contribution from the higher 
partial waves, since the numerical uncertainty of E^^\6d) 
is equal to 34 cm~^, and the numerical imcertainty of 
E^'^\7s) is equal to 1 cm^^. The numerical uncertainty 
of the second-order energy calculation for all other states 
ranges from 1 cm~^ to 5 cm~^. We use Z^ax = 6 in 
our third-order and all-order calculations, owing to the 
complexity of these calculations. Therefore, we use our 
high-precision calculation of £'^-' described above to ac- 
count for the contributions of the higher partial waves, 

i.e. we replace E^'^'> [^max = 6] with the final high-precision 

(2) (3) 
second-order value E^^^y. The contribution i?oxtra given 

in TableUaccounts for that part of the third-order MBPT 

correction not included in the SD energy. The values of 
(3) 

^cxtra ^'''6 quite large and including this term is impor- 
tant. 

We find that the correlation corrections to energies are 
especially large for 5/ states. For example, E^'^^ is about 
15% of and £(3) jg ^bout 36% of £(2) foj. 5J states. 



Despite the evident slow convergence of the perturbation 
theory expansion, the 5/ energy from the third-order 
MBPT calculation is within 0.9% of the measured en- 
ergy. The correlation corrections are so large for the 5/ 
states that inclusion of correlation leads to the different 
ordering of states using the DF and £ energies. If we 
consider only DF energies the 6^3/2 appears to be the 
ground state for Th IV, but the full correlation shows 
that the 5/5/2 is the ground state. The correlation cor- 
rections are much smaller for all other states; the ratios 
of £(0) and £(2) are equal to 6%, 5%, and 2% for the 6d, 
7s, and 10s states, respectively. 

The third-order and all-order results are compared 
with experimental values in Table |IT1 The energies are 
given relative to the ground state to facilitate compar- 
ison with experiment. Experimental energies for Fr I, 
Ac III, Th IV, and U VI are taken from and ener- 
gies of Fr-like Ra are taken from the NIST compilation 
[6[. Differences of our third-order and all-order calcu- 
lations with experimental data, 5E^^^ — E^^l — Ecxpt 
and SE^^ = E^^ — £cxpt, respectively, are given in the 
two final columns of Table [TTl In general, the SD results 
agree better with the experimental values than the third- 
order MBPT values. Exceptions are the cases where the 
third-order fortuitously give results that are close to ex- 
perimental values. Comparison of results from two last 
columns of Table |TT] shows that the ratio of SE'-^^ and 
SE^^ is about three for the 5/ states. As expected, 
including correlation to all orders led to significant im- 
provement of the results. Better agreement of all-order 
energies with experiment demonstrates the importance 
of the higher-order correlation contributions. 

Below, we describe a few numerical details of the cal- 
culation. We use the B-spline method described [13] to 
generate a complete set of basis DF wave functions for 
use in the evaluation of the MBPT expressions. We use 
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TABLE 11: The total third-order E^^^ and aU-order £fo? re- 
sults for Fr-like ions are compared with experimental energies 
Scxpt %&E = Etot - Sexpt. The energies are given relative 
to the ground state to facilitate comparison with experiment. 
Units: cm~^ 



nlj 


^ tot 


^ tot 


-fi'expt 


SE^'^ 


5E^° 






Fr I, Z = 87 








7Sl/2 

















7pi/2 


11683 


12156 


12237 


-554 


-81 


7P3/2 


13224 


13838 


13924 


-700 


-86 


6d3/2 


16282 


16048 


16230 


52 


-192 


6^5/2 


16251 


16217 


16430 


-179 


-213 


8Sl/2 


18893 


19695 


19733 


-840 


-38 


8pi/2 


22220 


23037 


23113 


-893 


-76 


8P3/2 


22737 


23575 


23658 


-921 


-83 


7d3/2 


23418 


24149 


24245 


-827 


-96 


7d5/2 


23461 


24238 


. 24333 


-872 


-95 




Ra II, Z ^ 88 








7Si/2 

















6^3/2 


12072 


12020 


12084 


-12 


-64 


6^5/2 


13460 


13561 


13743 


-283 


-182 


7Pl/2 


20936 


21182 


21351 


-415 


-169 


7P3/2 


25581 


25989 


26209 


-628 


-220 


5/5/2 


50182 


48728 


48988 


1194 


-260 


5/7/2 


50168 


48944 


49272 


896 


-328 




Ac III, Z = 89 








7Sl/2 

















6^3/2 


592 


825 


801 


-209 


24 


6^5/2 


3715 


4041 


4204 


-488 


-163 


5/5/2 


25701 


24018 


23454 


2247 


564 


5/7/2 


28132 


26420 


26080 


2052 


340 


7pi/2 


29110 


29303 


29466 


-356 


-163 


7p3/2 


37465 


37816 


38063 


-598 


-247 
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4138 
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4325 


-187 


-190 


6^3/2 


7640 


8304 


9193 


-1553 


-889 


6^5/2 


12646 


13377 


14486 


-1841 


-1109 


7Sl/2 


21891 


22229 


23131 


-1240 


-901 


7pi/2 


58732 


59213 


60239 


-1507 


-1026 


7p3/2 


71307 


71932 


73056 


-1749 


-1124 




U VI, Z = 92 








5/5/2 

















5/7/2 


7449 


7481 


7609 


-160 


-128 


6^3/2 


90359 


90593 


91000 


-641 


-407 


6^5/2 


99584 


99841 


100510 


-926 


-669 


7Sl/2 


141268 


141103 


141447 


-179 


-344 


7Pl/2 


192713 


191989 


193340 


-627 


-832 


7P3/2 


214669 


211747 


215886 


-1217 


-2689 
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FIG. 3: Coulomb second- and third-order energies as func- 
tions of Z for the 5fj and 6dj states in Fr-like ions. 
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FIG. 4: Valence removal energies {E/{Z — 84)'^ in cm ^) as 
functions of Z for the 5fj, 6dj, and 7si/2 states in Fr-like ions. 



negligible. 



50 splines of order fc = 8 for each angular momentum. 
The basis orbitals are constrained to a spherical cavity 
of radius R ~ 90-30 a.u for Fr I-U VI. The cavity radius 
is chosen large enough to accommodate all nlj orbitals 
considered here and small enough that 50 splines can ap- 
proximate inner-shell DP wave functions with good pre- 
cision. We use 40 out 50 basis orbitals for each partial 
wave in our third-order and all-order energy calculations, 
since contributions from the highest-energy orbitals are 



A. Z dependence of energies in Fr-like ions 

In Fig.m we illustrate the Z-dependence of the second 
and third-order energy corrections for the valence 
7s, 7p, 6d, and 5/ states of Fr-like ions. The second- 
order energy is a smooth function of Z for the 7pi/2 
and 5/7/2 states, but exhibits a few sharp features for 
the 7si/2 and 6^5/2 states. These very strong irregular- 
ities occur for the 7si/2 state Z=96, 99 and 6d^/2 state 
for Z=99 and are explained by accidentally small energy 
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denominators in the MBPT expressions for the correla- 
tion corrections to the energy. The third-order energy 
is a smooth function of Z for the 6d and 5/ states, 
but exhibits a similar sharp features for the 7pi/2 state. 
Most of the sharp features for the 6dj and 7pj states oc- 
cur at very high values of Z, Z > 97. Comparison of 
the and E'-^'> corrections for the 6d and 5/ states is 
illustrated by Fig. [31 This figure shows also the smooth 
dependence of second- and third-order corrections to the 
6^5/2 — 6(^3/2 and 5/7/2 — 5/5/2 fine-structure intervals. 

The total e'^^^^ energies divided by {Z - 84)^ in Fr-like 
ions are shown in Fig. [H We plot the five energy levels 
for the 7si/2, 6(ij, and 5/j states. The Z-dependence 

(3) 

of -E tot is smooth up to very high Z. Irregularities ob- 
served in Fig. [4] for high Z are explained by vanishing 
energy denominators in MBPT expressions for correla- 
tion corrections. 



III. ELECTRIC-DIPOLE MATRIX ELEMENTS, 
OSCILLATOR STRENGTHS, TRANSITION 
RATES, AND LIFETIMES IN FR-LIKE IONS 

A. Electric-dipole matrix elements 

The matrix element of a one-particle operator Z is 
given by [20| 



Zwv 



(5) 



where is the exact wave function for the many-body 
"no-pair" Hamiltonian H. In MBPT, we expand the 
many-electron wave function in powers of Vi as 



(0) 



(2) 



(6) 



The denominator in Eq. ([5]) arises from the normalization 
condition that contributes starting from third order [23] • 
In the lowest order, we find 



'(1) - 



(0) 



(0)\ _ 



(7) 



where z^v is the corresponding one-particle matrix ele- 
ment d^. Since "^w^ is a DF function we designate Z^^'> 
by Z(DF) below. 

The second-order Coulomb correction to the transition 
matrix element in the DF case with V^~^ potential is 
given by [l^ 



7(2) ^ ^ Zan{,Qwnva Qwnav^ ^ (.9' 



Jwavn ywanv l^na 



~t~ ^ ^ ~H Sn £y 

na na 

(8) 

Second-order Breit corrections are obtained from Eq. ([5]) 
by changing gijki to bijki, where bijki is the matrix ele- 
ment of the Breit operator given in (soj . 



In the all-order SD calculation, we substitute the all- 
order SD wave function ^'^'^ into the matrix element ex- 
pression given by Eq. ^ and obtain the expression [23| 



,(SD) _ ^wv + Z'-"-^ + ■ ■ ■ + 



^{1 + N^){1 + N,) 



(9) 



where Zwv is the lowest-order (DF) matrix element given 
by Eq. and the terms ZC"), k = a ■ ■ -t are linear 
or quadratic functions of the excitation coefhcients in- 
troduced in Eq. (jj]). The normalization terms Nw are 
quadratic functions of the excitation coefficients. As a 
result, certain sets of many-body perturbation theory 
terms are summed to all orders. In contrast to the energy, 
all-order SD matrix elements contain the entire third- 
order MBPT contribution. 

The calculation of the transition matrix elements pro- 
vide another test of the quality of atomic-structure cal- 
culations and another measure of the size of correlation 
corrections. Reduced electric-dipole matrix elements be- 
tween low-lying states of Fr-like systems with Z = 87-92 
calculated in various approximations are presented in Ta- 
ble |TTI1 

Our calculations of the reduced matrix elements in the 
lowest, second, and third orders are carried out follow- 
ing the method described above. The lowest order DF 
value is obtained from Eq. ([7]). The values Z'^^^^^^ are 
obtained as the sum of the second-order correlation cor- 
rection Z'^^^ given by Eq. ([5]) and the DF matrix ele- 
ments Z'^^^\ The second-order Breit corrections B^^'' 
are rather small in comparison with the second-order 
Coulomb correction Z^^^ (the ratio of B^-^^ to Z^^^ is 
about 0.2%~2%). 

The third-order matrix elements Z^-^^^^^^^ include the 
DF values, the second-order Z'^^^ results, and the third- 
order Z'^' correlation correction. Z'^' includes random- 
phase-approximation terms (RPA) iterated to all orders 

We find correlation corrections Z(^+^' to be very large, 
10-25%, for many cases. All results given in Table IIIII 
are obtained using length form of the matrix elements. 
Length-form and velocity-form matrix elements differ 
typically by 5-20% for the DF matrix elements and 2- 
5 % for the second-order matrix elements in these calcu- 
lations. 

Electric-dipole matrix elements evaluated in the all- 
order SD approximation are given in columns labeled 
Z(SD) (Eq. ^) of Table [ml The SD matrix elements 
^(SD) jjiclude Z^^^f completely, along with important 
fourth- and higher-order corrections. The fourth-order 
corrections omitted from the SD matrix elements were 
discussed recently by Derevianko and Emmons [3l| . The 
Z(SD) values are smaller than the Z(^^+^) values and 
larger than the Z^^^~^^~^^^ values for all transitions given 
in Table [ml 

In Fig. [51 we illustrate the Z-dependences of the line 
strengths for the 7si/2 — 7pj, 6dj — Ipj', and Qdj — 5fj' 
transitions. Two sets of line strengths values S'^^^ and 
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TABLE III: Reduced electric-dipole matrix elements in Pr-like ions calculated to first, second, third, and all orders of MBPT. 



Transition 


^(DF) 


^(DF+2) 


^(DF+2+3) 


^(SD) 


Transition 


^(DF) 


^(DF+2) 


^(DF+2+3) 


2-(SD) 






Pr I, Z = 87 










Th IV, Z = 90 








7Pl/2 


5.1438 


4.7412 


4.1362 


4.2641 


'51/2 


7Pl/2 


2.8994 


2.3748 


2.3669 


2.4196 


' '^1/2 


'P3/2 


7.0903 


6.6001 


5.6451 


5.8619 


( ^1/2 


'P3/2 


3.9933 


3.3399 


3.2930 


3.3677 




'Pl/2 


9.2216 


8.7781 


7.1945 


7.0276 


0^3/2 


'Pl/2 


2.5465 


2.1374 


2.0723 


2.1220 


0"3/2 


7p3/2 


4.2832 


4.1100 


3.3108 


3.2217 


0^3/2 


7p3/2 


0.9963 


0.8784 


0.8270 


0.8488 


oas/2 


'P3/2 


12.8041 


12.2858 


10.1204 


9.9456 


oa5/2 


'P3/2 


3.1975 


2.8348 


2.7006 


2.7549 


0"3/2 


5/5/2 


11.4529 


11.3002 


6.7352 


6.9611 


DCt3/2 


5/5/2 


2.4281 


1.4501 


1.3367 


1.5295 


AW 

0«5/2 


0/5/2 


3.0143 


2.9714 


1.9006 


1.9515 


0^5/2 


0/5/2 


0.6391 


0.4070 


0.3624 


0.4116 


0"5/2 


0/7/2 


13.4863 


13.2942 


8.5062 


8.7332 


005/2 


0/7/2 


2.9557 


1.8844 


1.7032 


1.9190 


Ra II, Z = 9,9, 






Pa V, Z = 91 






'Sl/2 


'Pl/2 


3.8766 


3.3758 


3.1793 


3.2545 


7Sl/2 


'Pl/2 


2.6267 


2.1111 


2.1365 


2.1822 


7Sl/2 


7p3/2 


5.3395 


4.7241 


4.3881 


4.5106 


7Sl/2 


'P3/2 


3.6173 


2.9720 


2.9750 


3.0327 


Oa3/2 


'Pl/2 


4.4462 


3.9314 


3.3975 


3.5659 


D03/2 


'Pl/2 


2.1785 


1.8079 


1.7817 


1.8152 


003/2 


7P3/2 


1.8815 


1.7070 


1.4262 


1.5117 


003/2 


'P3/2 


0.8319 


0.7327 


0.7004 


0.7119 


005/2 


'P3/2 


5.8616 


5.3340 


4.6240 


4.8232 


005/2 


'P3/2 


2.6895 


2.3817 


2.2936 


2.3185 


003/2 




5.3548 


4.6786 


4.5967 


4.4491 


003/2 


0/5/2 


1.9275 


1.0904 


1.0864 


1.2092 


DO5/2 


0/5/2 


1.4797 


1.3062 


1.3088 


1.2465 


005/2 


0/5/2 


0.5023 


0.3061 


0.2947 


0.3241 


005/ 2 


0/7/2 


6.6382 


5.8609 


5.8347 


5.6357 


005/2 


0/7/2 


2.3242 


1.4102 


1.3673 


1.5020 


Ac III, Z = 89 






U VI, Z = 92 






7si/2 


7pi/2 


3.2787 


2.7544 


2.6859 


2.7463 


7Sl/2 


7pi/2 


2.4165 


1.9138 


1.9562 


1.9599 


7Sl/2 


7p3/2 


4.5157 


3.8665 


3.7271 


3.8176 


7Sl/2 


7p3/2 


3.3271 


2.6952 


2.7168 


2.2449 


6rf3/2 


7pi/2 


3.1529 


2.6942 


2.5241 


2.6048 


6^3/2 


7pi/2 


1.9249 


1.5837 


1.5688 


1.5613 


6d3/2 


7p3/2 


1.2726 


1.1300 


1.0275 


1.0662 


6d3/2 


7p3/2 


0.7202 


0.6353 


0.6054 


0.4963 


6^5/2 


7p3/2 


4.0423 


3.6067 


3.3438 


3.4383 


6d5/2 


7p3/2 


2.3417 


2.0760 


1.9887 


1.6018 


6(^3/2 


5/b/2 


3.5764 


2.3879 


1.7460 


2.1624 


6^3/2 


5/5/2 


1.6274 


0.8798 


0.9034 


0.9974 


6^5/2 


5/b/2 


0.9577 


0.6708 


0.4760 


0.5877 


6^5/2 


5/5/2 


0.4209 


0.2479 


0.2470 


0.2674 


6^5/2 


5/7/2 


4.4148 


3.1076 


2.2956 


2.7539 


6^5/2 


5/7/2 


1.9492 


1.1353 


1.1293 


1.2293 



TABLE IV: Comparison of length [L] and velocity [V] results for reduced electric-dipole matrix elements in first, second, and 
third orders of perturbation theory in Fr-like systems with Z = 87-92. 



^(DF) ^{DF+2) ^(DF+2+3) ^(DF) ^(DF+2) ^(DF+2+3) 







L 




V 


L 




V 




L 




V 




L 


V 


L 




V 




L 




V 


87 


11. 


,4529 


11. 


,4002 


- 5/5/2 
11.2693 


11, 


.2693 


7, 


.8982 


7, 


,8972 


13, 


,4863 


6^5/2 

13.4022 


— 5/7/2 
13.2703 


13 


,2703 


9 


.6464 


9, 


,6455 


88 


5, 


.3549 


5, 


,1362 


4.7603 


4 


.7603 


4 


.4885 


4, 


,4839 


6, 


,6382 


6.3283 


5.9720 


5, 


,9720 


5 


.7071 


5, 


,7056 


89 


3, 


.5764 


3, 


,1632 


2.6134 


2, 


.6134 


0, 


.5254 


0, 


,4982 


4, 


,4148 


3.8730 


3.3832 


3 


,3831 


0, 


.8484 


0, 


,8576 


90 


2, 


.4281 


0, 


,7698 


1.6597 


1, 


.6595 


2, 


.8562 


2, 


,9000 


2, 


,9557 


0.5609 


2.1257 


2, 


,1255 


3 


.3454 


3, 


,3138 


91 


1, 


.9275 


2, 


.0349 


1.2702 


1, 


.2702 


1, 


.4017 


1, 


.4122 


2, 


.3242 


2.4373 


1.6097 


1, 


.6097 


1, 


.7219 


1, 


,7100 


92 


1. 


.6274 


1. 


.6199 


1.0312 


1, 


.0312 


1, 


.0559 


1, 


.0613 


1, 


.9492 


1.9305 


1.2978 


1, 


.2978 


1, 


.3024 


1, 


.2943 


87 


9, 


.2216 


9, 


6^3/2 

,9449 


- 7pi/2 

8.7906 


8, 


.7906 


6 


.8407 


6, 


.8365 


7, 


,0903 


7S3/2 

6.6425 


- 7p3/2 

6.6268 


6 


.6268 


5, 


.9488 


5, 


.9486 


88 


4, 


.4462 


2, 


,0213 


3.9815 


3 


.9816 


3 


.8547 


3, 


,8589 


5, 


,3395 


4.9435 


4.7968 


4, 


.7969 


4, 


.5150 


4, 


,5148 


89 


3, 


.1529 


2, 


,3683 


2.7533 


2 


.7533 


2 


.6221 


2, 


,6238 


4, 


,5158 


4.1785 


3.9641 


3 


.9641 


3 


.8018 


3, 


,8016 


90 


2, 


.5465 


2, 


,0531 


2.1960 


2 


.1961 


2, 


.1074 


2, 


,1084 


3, 


,9933 


3.7030 


3.4515 


3 


.4516 


3 


.3434 


3, 


,3432 


91 


2. 


.1785 


1. 


.8158 


1.8631 


1, 


.8631 


1, 


.7957 


1, 


,7964 


3, 


,6173 


3.3638 


3.0915 


3, 


.0915 


3 


.0121 


3, 


,0119 


92 


1. 


.9249 


1. 


.6384 


1.6343 


1, 


.6344 


1, 


.5743 


1, 


,5748 


3, 


,3271 


3.1029 


2.8190 


2, 


.8190 


2, 


.7460 


2, 


,7458 



5(1+2) presented for each transition. The values of 
5(1) and 5(1+2) obtained as [Z^^^^f and {Z^^'^+^^f, 
respectively. It should be noted that the values are scaled 
by {Z — 84)2 provide better presentation of the line 
strengths. The difference between 5(1) and 5(1+2) (.^py^y 
increases with increasing Z; for the 7s 1/2 — lp\/2 transi- 
tion, the ratio of the second-order (5(1+^) . 5(1)) and the 
first-order (5(i)) contributions is equal to 15% and 37% 



for Z =87 and 92, respectively. 



B. Form-independent third-order transition 
amplitudes 

Wc calculate electric-dipole reduced matrix elements 
using the form-independent third-order perturbation the- 
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TABLE V: Wavelengths A (A), weighted transition rates gA (s~^) and oscillator strengths gf in Ac III, Th IV, and U VI. 
The SD data (gyl'^^' and g/'^^') are compared with theoretical (gA*^™^' and g/'™'^') values given in Ref. m. Numbers in 
brackets represent powers of 10. 



Lower 


Upper 






g^{R.HF) 




^^(RHF) 








Ac III, Z = 89 








7Sl/2 


7p3/2 


2626.44 


1.63 [9] 


1.55[9] 


1.69[+0] 


1.6[+0] 


6^3/2 


7p3/2 


2682.90 


1.19 [8] 


1.23 [8] 


1.29[-1] 


1.3[-1] 


64/2 


7p3/2 


2952.55 


9.30[8] 


8. 23 [8] 


1.22[+0] 


1.1 [+0] 


7Sl/2 


7pi/2 


3392.78 


3.91[8] 


3. 59 [8] 


6.75[-l] 


6.2[-l] 


6^3/2 


7pi/2 


3487.59 


3.24[8] 


2. 75 [8] 


5.91[-1] 


5.0[-l] 


6^3/2 


5/5/2 


4413.09 


1.10 [8] 


2.34[8] 


3.22[-l] 


6.9[-l] 


6^5/2 


5/7/2 


4569.97 


1.61[8] 


3.02[8] 


5.04[-l] 


9.5[-l] 


64/2 


5/5/2 


5193.21 


4.99[6] 


1.03 [6] 


2.02[-2] 


4.2[-2] 








Th IV, Z = 90 








64/2 


7P3/2 


1565.86 


3.80[8] 


4.08[8] 


1.40[-1] 


1.5[-1] 


64/2 


7P3/2 


1707.37 


3.09[9] 


2.83[9] 


1.35[+0] 


1.2[+0] 


64/2 


7Pl/2 


1959.02 


1.21[9] 


1.04 [9] 


6.98[-l] 


6.0[-l] 


7Sl/2 


7P3/2 


2003.00 


2.86[9] 


2. 70 [9] 


1.72[+0] 


1.6[+0] 


7Sl/2 


7pi/2 


2694.81 


6.06[8] 


5.55[8] 


6.60[-l] 


6.0[-l] 


5/5/2 


64/2 


6903.05 


1.04[6] 


2.22[6] 


7.45 [-3] 


1.6 [-2] 


5/7/2 


64/2 


9841.54 


7.83[6] 


1.53 [7] 


1.14[-1] 


2.2[-l] 


5/5/2 


64/2 


10877.55 


3.68[6] 


7.93 [6] 


6.53[-2] 


1.4[-1] 








VYI,Z = 92 








64/2 


7p3/2 


800.729 


9. 72 [8] 


1.72 [9] 


9.34[-2] 


1.7[-1] 


64/2 


7p3/2 


866.737 


7.98[9] 


1.26[10] 


8.99 [-1] 


1.4[+0] 


64/2 


7pi/2 


977.129 


5.29[9] 


4.87[9] 


7.58[-l] 


6.9[-l] 


5/5/2 


64/2 


994.921 


1.47[8] 


3.51[8] 


2. 18 [-2] 


5.2[-2] 


5/7/2 


64/2 


1076.40 


2. 46 [9] 


5.55[9] 


4.26[-l] 


9.5[-l] 


5/5/2 


64/2 


1098.91 


1.52[9] 


3.64[9] 


2.75[-l] 


6.6[-l] 


7Sl/2 


7p3/2 


1343.39 


4.21[9] 


6.08[9] 


1.14[+0] 


1.7[+0] 


7Sl/2 


7pi/2 


1927.05 


1.09[9] 


1.03 [9] 


6.05[-l] 


5.8[-l] 



TABLE VI: Wavelengths A (A), weighted transition rates gA (s ^) and oscillator strengths gf in Fr I and Ra II. The SD data 
(5rA(^°) and 3/'^°') are compared with theoretical {gA^^™\ s/^^"^') values given in Ref. P, H and theoretical {gA^'^'^^'^\ 
gf(MBPT)^ values given in Ref. p3 ]. Numbers in brackets represent powers of 10. 



Lower 


Upper 






^^(MBPT) 


g^(RHF) 


ff/'^°^ 


^j(MBPT) 


^^(RHF) 










Fr I, Z = 87 










7Sl/2 


7P3/2 


7181.84 


1.88 [8] 






1.45[+0] 


1.46[+0] 


1.5[+0] 


7Sl/2 


7Pl/2 


8171.66 


6. 75 [7] 






6.76[-l] 


6.84[-l] 


6.5[-l] 


7pi/2 


74/2 


8328.18 


4.39[7] 






4.57[-l] 


4.40[-l] 


6.9[-l] 


7p3/2 


74/2 


9606.79 


8.83[7] 






1.22[+0] 


1.18[+0] 


l.l[+0] 


7p3/2 


74/2 


9689.14 


1.08 [7] 






1.52[-1] 


1.49[-1] 


1.2[-1] 


7Pl/2 


8si/2 


13342.0 


1.49[7] 






3.99[-l] 






7p3/2 


8Sl/2 


17216.1 


2.18[7] 


Ra II, Z = 88 




9.70[-l] 






64/2 


5/5/2 


2708.96 


2.02[9] 




1.59 [9] 


2.22[+0] 




1.7[+0] 


64/2 


5/7/2 


2813.76 


2.89[9] 




2.01 [9] 


3.43[+0] 




2.4[+0] 


64/2 


5/5/2 


2836.46 


1.38 [8] 




9. 80 [7] 


1.66[-1] 




1.2[-1] 


7Sl/2 


7P3/2 


3814.42 


7.42 [8] 


7.30 [8] 


7. 14 [8] 


1.62[+0] 


1.59[+0] 


1.5[+0] 


7Sl/2 


7Pl/2 


4682.28 


2.09[8] 


2.05[8] 


1.93 [8] 


6.87[-l] 


6.73[-l] 


6.3[-l] 


64/2 


7P3/2 


7077.95 


1.30 [7] 


1.29 [7] 


1.27[7] 


9.80[-2] 


9.70[-2] 


9.5[-2] 


64/2 


7P3/2 


8019.70 


9.13[7] 


9.10[7] 


7.85[7] 


8.81[-1] 


8.78[-l] 


7.6[-l] 


64/2 


7pi/2 


10788.2 


2. 05 [7] 


2.03[7] 


1.79 [7] 


3.58[-l] 


3.54[-l] 


3.1[-1] 
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FIG. 5: Line Strengths {S x {Z — 84)^ in a.u) as functions of Z in Fr-like ions. 



TABLE VIL Lifetimes r in ns for the nl levels for neutral fran- 
cium. Our SD results are compared with experimental mea- 
surements presented in Ref. fg| for the 7pj levels, in Ref. ^ 
for the 7dj levels, and in Ref. llj for the 8s 1/2 level. 



Level 


T-(SD) 


^(oxpt) 


7pi/2 


29.62 


29.45±0.11 


7P3/2 


21.28 


21.02±0.11 


7^3/2 


73.08 


73.6±0.3 


7^5/2 


67.93 


67.7±2.9 


8Si/2 


54.36 


53.30±0.44 



ory developed by Savukov and Johnson in Ref. [34] . Pre- 
viously, a good precision of this method has been demon- 
strated for alkali-metal atoms. In this method, form- 
dependent "bare" amplitudes are replaced with form- 
independent random-phase approximation ( "dressed" ) 



amplitudes to obtain form-independent third-order am- 
plitudes to some degree of accuracy. As in the case of 
the third-order energy calculation, a limited number of 
partial waves with /max < 7 is included. This restriction 
is not very important for considered here ions because 
third-order correction is quite small, but it gives rise to 
some loss of gauge invariance. The gauge independence 
serves as a check that no numerical problems occurred. 

Length and velocity- form matrix elements from DF, 
second-order (RPA), and third-order calculations are 
given in Table [TVl for the limited number of transitions in 
Fr-likc systems with Z = 87-100. The values differ 

in L and V forms by 2-15 % for the p — s transitions. The 
very large L — V difference (by a factor of 2-3) is observed 
in the Z'^^' values for d — f transitions as illustrated in 
Table IIVI The second-order RPA contribution removes 
this difference in L — V values, and the L and V columns 
with the Z(DF+2) headings are almost identical. There 
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are, however, small L — V differences (0.002%-0.2%) in 
the third-order matrix elements. These remaining small 
differences can be explained by limitation in the number 
of partial waves taken into account in the Imax in the 
third-order matrix element evaluations that we already 
discussed in the previous section describing the energy 
calculations. 



C. Oscillator strengths, transition rates, and 
lifetimes 

We calculate oscillator strengths and transition proba- 
bilities for the eight 7s — 7p, 7p — 6d, and 6d — 5/ electric- 
dipole transitions in Ra II, Ac III, Th IV, U VI and for 
the seven 7s — 7p, 7p — 7d, and 7p — 8s electric-dipole 
transitions in Fr I. Wavelengths A (A), weighted transi- 
tion rates gA (s~^), and oscillator strengths g f in Ac III, 
Th IV, U VI are given in Table|V]for Ac III, Th IV, U VI 
ions and in Table IVTl for Ra II and Fr I. 

The SD data {gA'^^^'> and are compared with 

theoretical (gA^"^^' and g/'™^)) results from Ref. Ij. 
The experimental energies were used to calculate the 
^^(SD)^ 5/(SD), gA(^^^\ and 5/(1^"") Therefore, we 
really compare the values of the electric-dipole matrix el- 
ements. The SD and RHF results for s — p and p — d 
transitions disagree by about 6-25 % with exception of 
the Qdj — 7p3/2 transitions where disagreement is 60 %. 
The largest disagreement (by a factor of 2-5) is observed 
between the SD and RHF results for the f — d transitions. 
The correlation corrections are especially large for these 
transitions as illustrated in Table [1111 Therefore, we ex- 
pect that our values, which include correlation correction 
in rather complete way, will disagree with RHF calcula- 
tions which appear not to include any correlation effects 
for transitions which involve 5/ states. Our conclusion 
is confirmed by comparison of the gA^^^^^ and gf'^^^^^ 
and our gA^^^^ and gf^^^^ results (see, also Ref p7}). 
Our values for transitions rates and oscillator strengths 
are in reasonable agreement (10-20 %) with RHF data. 

The SD data gA^^^^ and gf^^'^^ for Fr I and Ra II given 
Table IVTl are compared with theoretical data {gA^^'^^\ 
^y(RHF)) gi^gjj p,gfg^ j^jjj theoretical {gA'^^^^T^I 

and g/^^^^"^^) values from Ref. [ll|. The wavelengths 

;S^(0Xpt) g-^gjj Tg^l^lg [Y|] ^j^J^gjj JS^Jgrp 

pilation [6*] for Ra II and from Ref. [5] {E{7pi/2) = 
12237.409 cm-i, E{7p3/2) = 13923.998 cm-^), Ref. [| 
(£:(8si/2) = 19732.523 cm"!), and Ref. 5] (£^(7^3/2) 
= 24244.831 cm^^ and ^(7^5/2) = 24333.298 cm"!) for 
Fr I. Comparison of the gf and gA results obtained by 
three different approximations shows better agreement 
between our SD and MBPT [Ti] results than between our 
and HFR results obtained Biemont et al. in Refs. [H, [1], 
owing to more complete treatment of the correlation in 
our calculation and in Ref. [T^ . 

Our SD lifetime results are compared in Tablc lVlIl with 
experimental measurements presented in Refs. [8|, |9|, [Ol 
for the 7pj, 7c?j, 8si/2 levels in neutral francium. We 



find that our SD lifetimes are in excellent agreement with 
precise measurements provided in Refs. [1,11, [Il|. 



IV. MULTIPOLE MATRIX ELEMENTS, 
TRANSITION RATES, AND LIFETIMES IN 
FR-LIKE IONS 

Reduced matrix elements of the electric-quadrupole 
(E2) and magnetic-multipole (Ml, M2, and M3) oper- 
ators in lowest, second, third, and all orders of perturba- 
tion theory are given in Table IVini for Fr-like ions with Z 
— 88-92. Detailed descriptions of the calculations of the 
multipole matrix elements in lowest and second orders 
of perturbation theory were given in Refs. IH, f34| . 
Third-order and all-order calculations are carried out us- 
ing the same method as the calculations of El matrix 
elements. In Table IVlIIl we present E2, Ml, and M3 
matrix elements in the Z^^^\ Z^^^+'^\ ^(DF+2+3)^ ^^^^ 
Z^^'-*-' approximations for the &dj — 7si/2 transitions in 
Ra II and Ac HI and the Ml, E2 5/5/2 - 5/7/2 transi- 
tion in Th IV, Pa V, and U VI. We already mentioned 
that the ground state in Ra II and Ac III is the 7si/2 
state with the Qdj being the next lowest states; however, 
the ground state in Th IV, Pa V, and U VI is the 5/3/2 
state with the 5/5/2 being the next lowest state. The 
second-order contribution is about 1-3% for all transi- 
tions involving the 7si/2 states. For the 5/5/2 — 5/7/2 
transition, the second-order contribution (Coulomb and 
Breit) is very small (0.1%) for the Ml transition and is 
rather large (20%) for the E2 transition. 

In Fig. [HI we illustrate the Z-dependences of the line 
strengths for the 5/5/2 — 5/7/2, Qdj — 7si/2, and 5/j — 
7psi/2 transitions. Two sets of line strengths calculated 
in first- and second-order approximations are presented 
for each transition. The values of S^^^ and are 
obtained as {Z^^^^f and (Z(df+2))2^ respectively. The 
difference between S'^'^^ and 5'^-'^+^) curves increases with 
increasing Z . For the 5/5/2 — 5/7/2 transition, the ratio of 
the second-order (5(1+2) _ g{i)-^ ^^^^ first-order (S"'!') 
contributions is equal to 18% and 67% for Z =89 and 95, 
respectively. 

The strong irregularities occur in the curves describ- 
ing the second-order contributions (see, for example, 
5(1+2) (5_^^. _ 7^^^^)) fo^ z = and ^(1+2) (6^3/2 - 7si/2)) 
for Z — 95, 99). Those sharp features are explained by 
accidentally small energy denominators in MBPT expres- 
sions for correlation corrections as discussed above. 

Wavelengths and transition rates A'^^^^ for the elec- 
tric quadrupole (E2) and magnetic-multipole (Ml and 
M3) transitions in Ra II, Ac HI, Th IV, and U VI calcu- 
lated in the SD approximation arc presented in Table [IXl 
The largest contribution to the lifetime of the Qdj state 
in Ra II and Ac III ions comes from the E2 transition, 
but the largest contribution to the hfetime of the 5/7/2 
state in Th IV [l3| and U VI ions comes from the Ml 
transition. Our SD result for the Ml matrix element in 
Th IV \vf\ ion agrees to 0.5% with HFR results obtained 
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TABLE VIII: Reduced matrix elements of the electric-quadrupole (E2) and magnetic-multipole (Ml and M3) operators in first, 
second, third, and all orders of perturbation theory in Fr-like ions. 







Transition 


^(DF) 




^(DF+2) 


^{DF+2+3) 


^(SD) 








Ra II, Z = 


88 








M3 


6^5/2 


7Sl/2 


128.8000 




131.1100 


128.7695 


120.0000 


E2 


6^3/2 


7Sl/2 


17.2630 




17.0350 


13.7518 


14.5885 


E2 


6^5/2 


7Sl/2 


21.7710 




21.5580 


17.8089 


18.6906 


Ml 


6^3/2 


7Sl/2 


0.0000 




0.0008 


0.0189 


0.0014 




Ac III, Z = 


89 








M3 


6^5/2 


7Sl/2 


80.7820 




82.9520 


81.2050 


76.2511 


E2 


6^3/2 


7Sl/2 


10.6820 




10.4410 


9.2237 


9.5149 


E2 


6^5/2 


7Sl/2 


13.6550 




13.4270 


11.9668 


12.2808 


Ml 


6^3/2 


7Sl/2 


0.0000 




0.0015 


0.0203 


0.0013 




Th IV, Z = 


90 








Ml 


0/5/2 




1 R'^nc 

-L .OOUU 








1 A 


E2 


5/5/2 


5/7/2 


1.5669 




1.2339 


0.9724 


1.0834 






Pa V, ^ = 


91 








Ml 


5/5/2 


5/7/2 


1.8505 




1.8520 


1.8421 


1.8513 


E2 


5/5/2 


5/7/2 


1.1348 




0.8203 


0.7033 


0.7613 




U VI, Z = 


92 








Ml 


5/5/2 


5/7/2 


1.8505 




1.8518 


1.8443 


1.8512 


E2 


5/5/2 


5/7/2 


0.9153 




0.6205 


0.5413 


0.5864 



TABLE IX: Wavelengths A (A) and transition rates yl*^-^' 
(s^^) of the electric-quadrupole (E2) and magnetic-multipole 
(Ml and M3) transitions in Ra II, Ac III, Th IV, and U VI 
calculated in the SD approximation. Numbers in brackets 
represent powers of 10. 

Transition A A'^"' A 

Ra II, Z = 88 Ac III,Z = 89 

E2 7si/2 6d3/2 8275.145 1.536[-H00] 124844. 1.053[-06] 

E2 7si/2 6d5/2 7276.373 3.197[-H00] 23787.4 3.696[-03] 

M3 7si/2 6d5/2 7276.373 6.988[-07] 23787.4 7.070[-ll] 

Ml 7si/2 6d3/2 8275.145 2.281 [-05] 124844. 4.354[-13] 

Th IV, Z = 90 U VI, Z = 90 

Ml 5/5/2 5/7/2 23119.6 9.352[-01] 13143.0 5.089[-f00] 

E2 5/5/2 5/7/2 23119.6 2.487[-05] 13143.0 1.227[-04] 



by Biemont et al. in Ref. 1] since the correlation is small 
for Ml transition. 

Finally, we find that the lifetimes of the 6(^3/2 state is 
equal to 0.651 s in Ra II and 9.50x10^ s in Ac III; the 
lifetime of the 6(^3/2 state is equal to 0.313 s in Ra II and 
271 s in Ac III. The lifetime of the 5/7/2 state is equal to 
1.07 s in Th IV and 0.196 s in U VI. 



V. GROUND STATE STATIC 
POLARIZABILITIES FOR FR-LIKE IONS 

The static polarizability of Fr-like ions can be calcu- 
lated as the sum of the polarizability of the ionic core Ofc, 
a counter term a^c compensating for excitations from the 
core to the valence shell which violate the Pauli principle. 



and a valence electron contribution a^: 

a = ac + + Oi^c- (10) 

These contributions are given by formulas listed, for ex- 
ample, in Refs. [H,!!^. We calculate in the relativistic 
RPA approximation (see Ref. |37l] ) . 

The 7si/2 is the ground state in the cases of Fr I, Ra II, 
and Ac III, and the corresponding polarizability terms 
are given by 

N 

Qfv = ^[Iv{npii2) + Iv{np:i/2)], (11) 
n=7 
6 

n=2 

where 

I^(^1^)- .J^Uf"^''1 ' (12) 
3(2:; + 1) Enij ~ Ey 

and N is the size of B-spline basis set (N = 50 in this 
calculation). The calculation of the ay is divided into 
two parts: 

k 

^main ^ ^ [/, (npi/2) -f /. (71^3/2)] ; 
11=7 
N 

af^ = MnPi/2) + Iv{np3/2)]- (13) 

Here, k is equal to 10, 9, and 8 for Fr I, Ra II, and 
Ac III. The values of ajf^™ are calculated using SD val- 
ues of dipole matrix elements Zy^nij and experimental 
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FIG. 6: Multipole line strengths {S in a.u) as functions of Z in Fr-like ions. 



energies where they are available. We use SD energies 
when we did not find experimental data. The a*'^'^. and 
ttc contributions are small and are calculated in DF ap- 
proximation. 

Our numerical results are given in Table 1x1 The sum 
over n in the main polarizability term given by Eq. (|13p 
converges faster for Ac III than for Fr I. The ratios of the 
second and first terms in sum over n are equal to 1% for 
Fr I and only 0.1% for Ac III. Therefore, fast convergence 
allows us to limit the number of n in Eq. ([T5)) to n = 8 in 
Ac III, since we have no experimental energy values for 
high n for this ion. 

Our SD results for af^ given in last line of Table 1X1 are 
in good agreement with recommended value (317. 7±2. 4) 
for Fr given by Derevianko et al. in Ref. [38[ and rec- 
ommended value (104.0) for Ra II given recently by Lim 
and P. Schwerdtfeger in Ref. [3!|. We did not find any 
data for the a^^ in Ac III. Our ionic core polarizabili- 



ties ttc given in Table |X| are in an excellent agreement 
with recommended values (20.4 in Fr I and 13.7 in Ra II) 
presented in Refs. [38l.l39j. 

The valence polarizability for the 6/5/2 state, which is 
the ground state in the case of Fr-like ions with Z > 90, 
is given by [17] 

N 

"v = ^ [^("-^^3/2) + Ivind^/2) + Iv{ng7/2)]- (14) 

n—UQ 

Here, no equal to 6 for the nd states and 5 for the ng 
states. 

We use the same designations as we use for the 7si/2 
polarizability given by Eqs. pH]) - Our results are 

listed in Table |Xll The third-order (DF-h2-H3) and SD 
dipole matrix elements (a.u.) are calculated with the 50 
splines and cavity radius i? = 45 for Th IV. We use use 
experimental energies Q to calculate a™**™ in Table IXII 
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Both, third-order and all-order results for dipole polariz- 
ability of Th IV in the ground state 5/5 ji are presented 
in Table IXll (see also Ref. [l3|)- We see from this table 
that the largest contribution to the term comes 

from the Gd states (a™''™(6(i)). The core contribution 
(etc — 7.750 ao'^) is larger than the main term, limiting 
the accuracy of our calculations. The a*'''' and a^^ terms 
calculated in the DF approximation contribute only 5% 
to the final results. No experimental data are available 
for Th IV polarizability. 

The 5/5/2 ground state polarizabilities for Fr-likc ions 
with Z — 91-100 are given in Table IXlTl Results are ob- 
tained in DF and second-order MBPT approximations. 
The contribution of the 537/2 state into the as/^^^ polar- 
izability increases with increasing Z\ 1% for ion with Z 
= 90 (Th IV), 19% for ion with Z = 92 (U VI), and 30% 
for ion with Z = 95 (Am IX). The main contributions 
are smaller than the core contributions. The calcula- 
tions are conducted only in second-order approximation 
for these ions, owing to problems with accidentally very 
small denominators in the corresponding calculations of 
the properties with involving hg state. We mentioned 
previously that the core values ac(c) are calculated in 
the relativistic RPA approximation following by method 
described by Johnson et al. in Ref. js^]. Our ac(c) re- 
sults for Ra II - U VI disagree with results presented by 
Biemont et al. in Ref. [l|, but agree with Ra II result 
given in Ref. [3§|. 



VI. CONCLUSION 

In summary, a systematic relativistic MBPT study of 
the atomic properties of the 7si/2, 7pj, Qdj, and 5/j 
states in Fr-like ions with nuclear charges Z = 87 — 100 
is presented. The energy values are in good agreement 
with available experimental energy data and provide a 
theoretical reference database for the line identification. 
A systematic all-order SD study of the reduced matrix 
elements and transition rates for eight 7s — 7p, 7p — 6(i, 
and 6d — 5f electric-dipole transitions is conducted. Mul- 
tipole matrix elements (7si/2 — Sdj, 7si/2 — 5/j, and 
5/5/2 ~ 5/7/2) are evaluated to obtain the lifetime data 
for the 6(i3/2 and the 5/7/2 excited state. The scalar po- 
larizabilities for the 7si/2 ground state in Fr I, Ra II, and 
Ac III and 5/5 /2 ground state in Th IV are calculated us- 
ing a relativistic third-order and all-order methods. The 
scalar polarizabilities for Fr-like ions with nuclear charge 
Z = 90-100 in the 5/5/2 ground state are calculated us- 



ing a relativistic second-order MBPT. These calculations 
provide a theoretical benchmark for comparison with ex- 
periment and theory 
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TABLE X: Contribution to the 7si/2 static polarizabilities 
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Fr I 


Ra II 


Ac III 
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46.74 
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0.23 


0.06 




0.54 


0.03 




a--"(10p) 


0.19 






main 
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„ tail 


1.20 


0.11 


0.04 
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20.40 


13.79 


11.42 


Give 


-0.95 


-0.74 


-0.54 




313.7 


106.5 


57.71 



TABLE XL Contribution to the 5/5/2 static polarizabilities 
(a.u.) of Th IV calculated in the third-order - (a) and all- 
order - (b) approximations. 



v(«) 



"(6d) 
"(7d) 
"(8d) 
"(557/2) 

"(6<?7/2) 
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